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Abstract 

The Vlasov-Poisson system describes interacting systems of collision- 
less particles. For solutions with small initial data in three dimensions 
it is known that the spatial density of particles decays like t~ 3 at late 
times. In this paper this statement is refined to show that each derivative 
of the density which is taken leads to an extra power of decay so that in 
N dimensions for N > 3 the derivative of the density of order k decays 
like t~ N ~ k . An asymptotic formula for the solution at late times is also 
obtained. 



1 INTRODUCTION 



The Vlasov-Poisson system provides a statistical description of the dynamics of a 
large number of particles which are acted on by a force field which they generate 
collectively. One class of applications of this system is in plasma physics where 
the force is electrostatic and the particles are electrons or ions [10]. Another is 
in stellar dynamics where stars play the role of particles. The particle treatment 
is justified in models of galaxies where the distance between stars is much larger 
than their diameters. In this case the force is gravitational [2]. The equations 
in these two cases only differ by a sign and a lot of the mathematical theory 
works in exactly the same way for both. This applies in particular to the results 
of this paper. For surveys of results on the Vlasov-Poisson and related systems 
see [8] and [1]. 

The distribution function / of the particles satisfies the Vlasov equation 
while the potential cf> for the field satisfies the Poisson equation. The function 
/ depends on time t, the spatial point x G M 3 and the velocity v € K 3 . It is 
natural to pose an initial value problem with / being prescribed at t = 0. For an 
initial datum which is C 1 and has compact support it is known that there exists 
a unique corresponding C 1 solution, globally in time [18], [15]. The support of 
/ is compact at each fixed time t and an important diagnostic quantity is P(t), 
the supremum of \v\ over the support of / at time t. Estimates are known for 
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P(t) [12] and these imply estimates for ||p(t)||j,oo where p, the spatial density of 
particles, is given by p(t, x) = J f(t, x, v)dv. We have P(t) < C(l +t) log(2 + i). 
Unfortunately these estimates seem far from optimal. They are the same for 
the plasma physics and stellar dynamics cases. Intuitively it is to be expected 
that the optimal estimates differ in these two cases. In the stellar dynamics 
case there exist time-independent solutions so that ||p(i)||z,°° does not decay in 
general. In the plasma physics case decay estimates for integral norms of p are 
established in [13] and [17]. The pointwise estimates can also be improved to 
give a bound for P(t) of the form C(l + tf/ z [21]. 

It is possible to consider the analogue of the Vlasov-Poisson system in higher 
dimensions. It is, however, known that global existence fails in four space di- 
mensions [11]. An explicit example of singularity formation and information on 
the asymptotics of solutions near a singularity were obtained in [14]. 

There is a case where much more is known about the long-time asymptotics 
of solutions of the Vlasov-Poisson system, namely that of small initial data. 
The first global existence theorem for that case due to Bardos and Degond [3] 
naturally comes with decay estimates. They show that 

<C(l + t)- 3 

If the data are sufficiently differentiable then the same techniques should lead 
to estimates of the form 

\\D k P (t)\\ L ~ <c(i+<r 3 

but they do not give more. In this paper we apply new techniques to this 
problem to obtain estimates of the form 

\\D k p{t)\\ L oo <C(l+ty 3 - k 

for solutions with small initial data. Furthermore, we obtain asymptotic expan- 
sions for these solutions. 

Note that there are a number of generalizations of the results of [3] in the 
literature. The fully relativistic generalization of the plasma physics problem is 
given by the Vlasov-Maxwell system. An analogue of the result of [3] in that 
case was proved in [9]. In the stellar dynamics problem the fully relativistic 
generalization is the Einstein- Vlasov system [24] which is much more compli- 
cated. A small data global existence theorem in the spherically symmetric case 
was obtained in [22]. A related system which is physically incorrect but mathe- 
matically interesting is the Vlasov-Nordstrom system for which there is a global 
existence theorem [5]. Surprisingly it seems that no analogue of the asymp- 
totic result of [3] has been proved for this system. There are generalizations 
of the results for solutions of the Vlasov-Poisson and Vlasov-Maxwell systems 
with small data to almost spherically symmetric data [25], [19]. There are also 
results for solutions of the Vlasov-Poisson system with non-standard boundary 
conditions which are relevant to cosmology [23], [20]. Global existence has also 
been proved for some cosmological solutions of the Einstein- Vlasov system with 
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symmetry. See for instance [26]. It would be interesting to extend the results 
of this paper to some of the cases mentioned in this paragraph. 

This paper was motivated by the wish to prove a small data global existence 
theorem for the Einstein- Vlasov system which does not require any symmetry as- 
sumptions. To understand the difficulty of this problem note first that even the 
vacuum Einstein equations, from the present point of view the Einstein- Vlasov 
system with / = 0, are very hard to handle mathematically. The landmark 
work of Christodoulou and Klainerman on small data global existence for the 
vacuum Einstein equations [6] is so complicated as to discourage any attempts 
to incorporate matter. The more recent alternative proof of Lindblad and Rod- 
nianski [16] looks much more promising. Nevertheless, it seems to require good 
decay estimates for higher derivatives, i.e. estimates similar to those proved for 
the Vlasov-Poisson system here. 

The Vlasov-Poisson system in N dimensions reads 

/ t + wVs/ + 7 V^-V„/ = , xtR N ,t>0 (1.1) 

A<f)= I fdv = p(x,t) , xeR n ,t>0 (1.2) 
Jr n 

where / = f(x,v,t). In the following we assume that f(x,v,0) = fo(x,v) 
has finite L 1 norm and N > 3. The sign 7 = ±1 corresponds to the plasma 
physics and gravitational problem respectively. Since the results of this paper 
apply equally to both cases, we will restrict our analysis to the case 7 = 1. No 
sign condition on / is needed. However, some additional decay properties for 
fo (x, v) will be assumed. 

Global existence and decay estimates for the Vlasov-Poisson system were 
studied in [3] in three spatial dimensions under suitable smallness and regularity 
assumptions for /q. These estimates are optimal in the rate of decay for the 
density p since, for small compactly supported initial data, the volume of the 
support of p can be bounded by C (1 + 1) 3 so that if the decay in L°° was 
stronger than (1 + i) 3 , the total number of particles (i.e. the L 1 norm of p) 
would decay, leading to a contradiction due to the conservation of that quantity. 
However, they do not provide the optimal rate of decay for the derivatives that 
could be expected on dimensional grounds. 

For small initial data the dynamics of the Vlasov-Poisson system might be 
expected to be dominated by the free streaming part of the equation: 

ft+v- V x / = 

because the term V x ■ V„/ is quadratic in the density. (Actually this is a 
consequence of the Bardos-Degond analysis) . If we assume that the dynamics 
of the problem is dominated by the free streaming regime as t — > 00 and the 
initial density of particles is, say, compactly supported (fast enough decay works 
similarly), the velocities of the particles would be bounded by a number of order 
one. Therefore, the support of the density p would spread linearly The field V(f> 
generated by a particle density with finite mass spread over a region of order t 
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decreases as ^ as can be easily seen by means of a rescaling argument. Notice 
that a posteriori this provides a justification for the assumption that was made 
before concerning the finitcness of the deviation of the velocities of the particles 
due to the interaction of the field. 

The main contribution of this paper is the development of a technique that 
allows us to obtain optimal decay estimates for the solutions of the VP system 
in iV-dimensional space. More precisely, the rescaling argument sketched above 
suggests that the particles spread into a region of volume t N in the x-coordinate. 
Since the total mass of the particles is of order one it would be natural to expect 
the following estimates for the density 



\p\< 



|Vp|< 



v P\ < 



c 



(*+l) 

c 



N 



(*+l) 

c 



JV+1 



(* + l) 



Af+2 



|V fe p| < 



c 



(*+l) 



N+k 



The first estimate was obtained by Bardos-Dcgond for the case N = 3 and 
can be similarly extended to the case N > 3. Our method allows us to obtain 
the corresponding estimates for the derivatives for small initial data. 

The basic idea of the method is as follows. It is easy to see self-similar be- 
haviour for the density (and the derivatives) in the free streaming case. Indeed, 
in that case, integration along characteristics yields 



/ {x, v, t) = f (x - vt, v) 



whence: 



p(x,t) = J f{x,v,t) dv = J f (x - vt, v) dv 
In order to obtain self-similar behaviour we make the change of variables 



whence: 



x = 
dv = 

p(x,t) 



o — vt 

d0t (&r ) 



dx = -jjfdx 



1 

w 



fo xo 



x 



t 



dxo 



In the limit t — > oo this formula yields the self-similar behaviour in the region 
where \x\ is of order t: 



p(x,t) 



J fo (xq, I) dx = -j^-Pfs (|) (1.3) 
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Here the asymptotic free streaming density pf s is given by 

PU (y) = J fo(xo,y)dx 

Notice that (1.3), at least formally, provides the desired estimates for the 
derivatives of p. The key idea of our argument is a method for generalizing this 
method to the full VP system with small initial data. The main point is the 
following. Suppose that the characteristics starting at xq, vq reach the points 
x, v at time t. Assuming suitable invertibility conditions any pair of variables 
in the set (xo, vq, x, v) can be used as a set of independent variables in order 
to represent the others. The previous argument for the free streaming case 
suggests using x, xo as independent variables. However, in order to determine 
the functions that provide Vo, v in terms of x, xq it turns out to be necessary 
to solve a boundary value problem for the characteristic equations. The main 
argument of this paper consists in proving that such a boundary problem can 
be solved for small initial densities and that the corresponding solutions of such 
a boundary value problem satisfy suitable regularity and decay estimates. 

Using a similar method it is possible to obtain not only estimates for the 
derivatives of the density, but also convergence of the solutions of the VP system 
to a self-similar solution. More precisely, we rewrite the problem using the self- 
similar variables y — ^pjy, v — v, r = log(t+ 1) , / = ^jyrc g and after 
integrating the resulting equations along characteristics we replace the variables 
(y, v) by (y, yo) . This change of variables requires the solution of a boundary 
value problem for the characteristic equations analogous to the one described 
above. Such a boundary value problem can be analyzed in detail as t — > oo, 
and this provides the asymptotic behaviour of the density p as t — > oo. One of 
the relevant results of the analysis is the fact that although the asymptotics of 
the solutions is self-similar, the precise function describing the asymptotics of 
the density depends in a very sensitive manner on the choice of the initial data 
fo (x, v) . This analysis is done in the last section of the paper. 

The paper is organized as follows. In Section 2, we derive estimates for the 
density and its derivatives. In Section 3, we prove convergence to the self-similar 
solution. Throughout the paper, C > will denote a generic constant that may 
change from line to line and is independent of t, e 0} fo- 

2 ESTIMATING p AND ITS DERIVATIVES. 

2.1 The main result. 

We will use the following function spaces extensively. 
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= sup{ / \p (x, t)\ d N x + (t + if V (i + 1) £ || VV (■, t) 

+ (t + l) 7 
where < a < 1. 



1=0 

,N+k +a \V k p(x,t)-X7 k p(y,t)\ 

SUP : ^ } 

^yGR" F _ V\ 



(2.1) 



fe+2 



\\4>\\ Yh a = sup{(t + if" 2 £ (t + 1) £ || VV (■, *)|U (KW) 

|V fe+2 </>(*,*)- V fe+ V(z',i)| 



(2.2) 



+ (*+!) 



Af-2+fc+a 



sup 



X — X 



} 



The main result of the paper is the following Theorem: 
Theorem 1 Suppose that f (x, v) satisfies the following assumptions: 



k I 

EE 



e=o m=0 



d e fi 







E 



sup 

~ n x,x'eR' 



Q x mg v £-rn 
dx m dv k ~ m ^ dx m dv k ~ m ^ ' ^ 



< 



(l + \x\) K (l + \v\f 



(2.3) 



< 



m=0 



\X — T 



(i + H) 



K 



< a < 1 



E sup 



a*7o 



dx m dv k - 



(x,v) 



d k fa 
dx m dv k ~ 



(x,V>) 



< 



£o 



\v — V 



(l + \x\f(l + \v\f 



, < a < 1 , 



for some suitable K > N and e > small enough. Then there exists a corre- 
sponding solution of the Vlasov-Poisson system with 



\x h , c 



< Ce 



A result analogous to Theorem 1 was proved in the case k = under slightly 
different assumptions on / , by Bardos-Degond (cf. [3]). The main contribution 
of this paper is to derive the optimal decay estimates for the derivatives of p. 
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2.2 A basic boundary value problem for the characteristic 
curves. 



We introduce some basic notation. Suppose that the characteristics starting at 
(xo,Vo) at time t = reach the point (x,v) at time t. The basic idea in the 
paper is to use x, x as independent variables to describe the values of v and 
vq. More precisely, we will write: 



V = W (x,Xq, t) 
V = Wq (x,X ,t) 



(2.4) 
(2.5) 



Note that the existence of the functions w, wo is assured by the implicit function 
theorem as well as the estimates that we will derive later where it will be shown 
that the change of variables is close to a change of variables that can be inverted 
explicitly. By changing the variable v to x in the integral with x and t fixed, 
it then follows, using (2.4) that 



dv 



det 



dw (x, xq, t) 



dx 



dxo 



whence 
p{x,t) = 



R N 



f(x,v,t) = / f(x ,v )dv 



J f(xo, 



Wq (x,XQ,t)) 



det 



dw (x, xo,t) 



dx 



We now formulate the following auxiliary boundary value problem that de- 
scribes the evolution of the characteristics starting at the spatial point xq at the 
initial time and reaching the point x at time t : 



"' ' ' ;; ) =V(a), = V<£ (X (a) , s) , X (t) = x, X (0) = x 



ds 



ds 



(2.6) 



Notice that the functions X (s) , V (s) depend also on the variables x, xq, t. 
However, for simplicity, we will not write the dependence on these variables 
explicitly unless it is needed. We then rewrite the above characteristics as a 
perturbation from those associated to the free streaming case as follows 



dX(s) 
ds 



V(s) 



X — Xq 
t 



dip 
ds 



V<j>(X{s),s), X(t)=x, X{0)=x , 

(2.7) 

where <p(s) — ip(s;x,xo,t) is the perturbed value of the velocity with respect 
to the free streaming case. Notice that in the limit of zero density p = 0, the 
field 4> vanishes and ip (s) = 0. 

We examine the derivatives of the function ip (s) with respect to the variables 
x, xq, in order to derive suitable estimates for p. The density function p can be 
represented as 



p(x,t) = J f(x,v,t)dv = J f {x ,V(0;x,w(t,x,x 



),*)) 



/ dw 
\dx 



dxo. 
(2.8) 



dxo. 
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Along the characteristics, we have 



dw dV , , 1 dip , , 

ax cteo £ cteo 

where Jjv is the A^-dimensional identity matrix. 

On the other hand, we wish to obtain estimates for the derivatives of p. 
Suppose for the moment that we restrict our attention to the first derivative of 
p with respect to x. Such a derivative is given by: 



dp 
dx 



(x,t) 



d f dV 
-^(x o ,V(0;x,w(t,x,x ),t)) — (0) 



( dw\ 
\dx a J 



dxo 



+ 



f d 

/ fo (xo, V (0; x, w (t, x, x ) , t)) — 



/ dw 
\dx 



dxo- 



To estimate the first derivative of p reduces to derive estimates for: 

d 2 v 



dV , dV . 



dxdxo 



(s = t). 



Equivalently 



Notice that the equation of the characteristics (2.7) indicates that in order 
to obtain bounds for two derivatives with respect to x, xq of the characteristic 
curves we need to estimate three derivatives of the potential </>. These are the 
exact number of derivatives that can be expected to be estimated from the 
Poisson equation under the assumption that |^ is bounded. Nevertheless, in 
order to avoid the standard problems that arise in the regularity estimates for 
the Poisson equation in the spaces C k , it is necessary to work with the Holder 
spaces C k ' a . 



2.3 Estimates on the regularity and the rate of decay of 
f (s; x, x , t) in terms of the properties of the potential 

We present a key a priori estimate for ip in terms of 4> m the following Propo- 
sition. We define two norms with respect to the spatial variable x. 

Definition 2 For u (•) E L°° (R N ) , 

II u IIl°o (x) = ess SU P \u(x)\, 



For u (•) G C°' Q (R N ) , 



sup 



- (xi) - u (x 2 )| 
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For notational simplicity, in the following, we will use \\u (s) H^oo^) , [u (s)] a 
instead of ||w (s; •, xq, t) || L oo/ x ) , [u (s; •, xq, t)] a ,y which in fact depend ons,so, 

and t. For example, ^f- (s) will denote 



: a,(x) \- UJj J0,a,(x) 

Proposition 3 Suppose that 

for a suitable e > sufficiently small. Suppose that t>l. Then, the following 
a priori estimate holds: 



t sup 

0<s<t 



dip 



+ t 



1+a 



f 

Jo 



L°°{x) 

1 dip 



+ t' 



0,a,(x) 



dx 

d 2 v 



(s) 



L°°(x) 



+ t 1+a sup 

0<s<t 



dip 
dx 



(s) 



<C\ 



0,a,(x) 



dxdxo 



(t) 



+ t 



2+a 



dx 



(s) 



ds + t 



+ t ds 
Jo 



L°o(x) 

d 2 <p 



7' 

Jo 



dp 
dx 



(s) 



L°°(x) 

ds 



d 2 v 

dxdxo 



(*) 



<C\ 



0,a,(x) 



^1, 



+ t 



1+a 



L°°(x) 



0,a,(x) 

d 2 y 

dxdx n 



(s) 



ds 



0,a,(x) 



dxdxo 

for some suitable constant C > independent oft, go- 
Proposition 3 is the main new technical result of the paper. This estimate 

provides optimal decay properties for the derivatives of the characteristics in 

terms of the decay properties of the derivatives of the potential 4>. 

Furthermore, we derive a generalization of Proposition 3 under additional 

regularity and decay assumptions for the potential <f> ■ 

Proposition 4 Suppose that 

H\\ Ytia < eo, £>2 

for a suitable e > sufficiently small. Suppose that t > 1. Then, the following 
estimates hold: 



Y,t k sup 



fc=i 



0<s<t 



d k p> 
dx^ {s) 



+ t e+a sup 



L°°(x) 



0<s<t 



0,a,(x) 



<cu\\ Ytti 



^ L 



t Qk+l 



dx k dx a 



ds + t e+a 



L°°(x) 



f 

Jo 



d e+1 p 

dx e dx 



(s) 



0,a,(x) 



k=i 



k+1 



Qk+l 



dx k dx 



(s) 



+ t 



t+l+a 



L^(x) 



Qt+1 



dx e dx a 



(*) 



(2.9) 

<is<cu\\ Yea 

(2.10) 

<C|Wy (i> 

(2.11) 



0,a,(x) 
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for some constant C > independent oft, Eq- 



2.3.1 Preliminary results: Integral equation satisfied by ip (s, x, Xo, t) . 

The perturbed velocity tp (s) satisfies the following integral equation 
Lemma 5 ip(s) = p (s;x,x ,t) in (2.7) satisfies the integral equation 



where 



G(o = V0(x(o,O. 



X{H)=x + X —^ti + j\(s) 



ds. 



(2.12) 

(2.13) 
(2.14) 



Proof. We integrate (2.7) to obtain 



<P(s) = <p(0) 



G(Ode, X(s) = x 



x - x 



The boundary condition X (t) = x yields 



/ p (s) ds = 0. 
Jo 



s+ / p(s)ds. (2.15) 



(2.16) 



Therefore, integrating the first equation in (2.15) from s = to s = t and using 
(2.16), we have 

1 



t 



ds. 



<p(0) 

Thus, we obtain from (2.15) 

p( S ) = -\f Q [f^GiOK 

By changing the order of integration in the first integral on the right-hand side 
of this formula, we deduce (2.12). ■ 



ds 



G (Ode- 



Note that 



V (t) = V (t; x, x , t) = v = w (x, x , t) = 

V (0) = V (0; x, x a , t)=v = w (x, x , t) 



x - x 

t 

X — Xq 

~ t 



p(t;x,x ,t) , 
<p{0;x,x o ,t) . 



The integral equation (2.12), (2.13), (2.14) is the key ingredient that will 
be used to derive optimal regularity and decay estimates for the functions 
w (x, x , t) , Vo (x, x , t) . As a first step we prove that the solutions of (2.12), 
(2.13), (2.14) are well defined if eo is small enough. 
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Lemma 6 (Solvability) Let k > be an integer. There exists Sq > such that, 
for any t > 1, x G R w , xo € R^ and any function <f> satisfying 



< £o, 



i/iere exists a unique solution ip(-) = <p (■; x, x ,t) e C([0,i]) of (2.12), (2.13), 
(2.14). 

Proof. Let the space of functions for ip (s) be 



X = ^C([0,t]) : sup |^(s)|<l . 

0<s<t 



Let 



where 



J M (*) = - jf G v (0 ^ + 7 £ ^ (0 ^. ( 2 - 17 ) 



(0 = x + 



X — Xq 

T 



We first show that J' is a well-dehned operator in the space X. By definition, 
we have 



|G V (OI<||V0(-,OII 



< 



This yields 



\J(<p)(s)\ < / 



< 



< 



\Y k „ 



\Y k „ 



(1 + 



JV-1 



1 + 



log (1 + t) 



tJo (1 + 



where we have estimated the last integral term using the fact that AT > 3. This 
idea will be used repeatedly in the following. Thus J (ip) is bounded for all 
t and J is well-defined in the space X. We next show that the operator J is 
contractive. Taking the difference of J (tf>\) and J (tpi) yields 



[J (<Pi)-J M (s) = -J* [G V1 (0 - G V2 (0] de+7 £ € [G V1 (0 - G V2 (£)] 



ii 



Using the definition of 



we have 



Sup \[J{<Pl)-J(<P2)](8) 
0<s<t 



< f \\V 2 <t>U)\\ L ~ ix) f \Ms)-Ms)\d8 dt + jj\\\V 2 t(;0\\ L ^ x) J*\<pi(a)-<p2(8)\d8 



< 



in.o 



Js (4 + 1 ; Jo 1 Jo + Jo 



(£+l) J 



< Ce sup |(yi - y 2 ) + C£ 1OS ^, + — sup |(<^i-y 2 )(s)| 

0<s<t t 0<s<t 

< Ce sup |(^i - ip 2 ) (s)| . 

0<s<t 

Thus we can choose £0 small enough such that Ceo < 1 in the above inequality 
and conclude that J is contractive. Notice that C is independent of t, x n ,x and 
eo can be chosen independently of these variables. Then by the Banach fixed 
point theorem, we deduce the existence and uniqueness for <p satisfying (2.12) 
in the space X. m 

Corollary 7 Let ip (s) be the solution in Lemma 6. Then we have 



f 

Jo 



<p (s) ds = 0. 



Proof. Using (2.17) we derive the identity: 



Jo 



J (</?) ds = 0. 



Therefore, using Lemma 6, we have 

<P = J(<P)- 

This completes the proof of Corollary. ■ 

Notice that this Corollary implies that the characteristics (2.7) satisfy the 
desired boundary condition for X (s) , i.e., X (t) = x. 

We now turn to decay estimates for the density function p and its derivatives. 
Before we proceed, we state some basic properties of the Holder norms. 

Lemma 8 

[/s] ,a,(x) ^ ^{ll/llioo [ff]o,a,(x) + [/]o,a,(x)} . 

for any f,g€ L°° nC°<", 

[/]o,a,(x)<C||/||i-«||V/||^, 
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for any f G TT 1,00 , 

F°4>,«,(s) <[*!(>,<,,(*) l|V<oo, 

/or any F G C 0,a and any u G jy 1 ' 00 . 

Proof. The results in the lemma are standard estimates for Holder norms. 



2.3.2 The proof of Proposition 3. 



The proof of Proposition 3 follows from a sequence of lemmas. There are three 
ideas that will appear repeatedly in all the remaining arguments of this paper. 
Estimating terms like -g^- and its derivatives with respect to x, it is not possible 
to obtain bounds for the rate of decay suggested by dimensional considerations 
for all the values of s G [0, t] . It is, however, possible to obtain such optimal 
decay estimates for the integrals of such terms in the interval [0, t] as well as 
for the time s — t which is the only one where such optimal estimates are 
really needed. The second idea is that it is convenient to obtain, before deriving 
pointwise estimates, integral estimates for terms like J Q (s) ds. The 

third idea is that the estimates for terms that do not contain derivatives with 
respect xq are more easily obtained by directly estimating the supremum over 
the interval [0, t] and using Gronwall-type arguments, without any need for 
estimating integrals over the interval [0, t]. 



Lemma 9 There exists e small such that for t > 1 and any function 4> satis- 
fying 

< eo, 



we have 



f 

Jo 





ds<cu\\ Yoa , 

L°°(x) 


(2.18) 


£«> 


Loo f x ) t 


(2.19) 



(2.20) 



Proof. Differentiating (2.12) with respect to xq yields 
where, for simplicity G v — G. We now take of (2.13)-(2.14) to get 

Lg (0 = W (x (0 , ^-x (0 = V* (x (o , o [(i + 



_d_ 

dxr. 
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and since 4 < 1, we have 



< 



l»w (£ + 1) 
Ml* 



JV 



< 



" K + i) 

Putting the above into (2.20) yields 



AT 



Jo 



9x 



L°°{x) 



< C U\\y , a 



K + l) 



JV 



[1 



dx 



(2.21) 



Z,oo (x) 



+ 



ox 



L°°(x) 



By integrating the above from to t and by the assumption, we obtain 



dip 
dx 



to 



ds<CU\\ Yoa f [*—*[!+ f 



dp 



to 



+ 



' Jo 



<cu\\ Y +cu\\ Yo 



N 



dx 



o 



9x 



ds 



L°°(x) 

ds}d£ 



ds]d£ 



<C\ 



+ C*e 



<9</5 



to 



L°°(x) 

ds, 



where we have used the estimate J Q 



and the fact that J Q * 



ds = Jo 

small enough so that Ce < 1/2, we get 



t d£ 



* 1 



L~(x) 



^</o J^to 



Z,~(x) 



(is 



Jo ^ ^ < C. Thus if e is 



/' 

Jo 



|^to 



ds<CU\U 



(2.22) 



L°°(x) 



Putting s = f in (2.21) and using (2.22) yields 



L°°(x) 



to 

OT 



L°°(x) 



<c^(i + w y j<c . 



Y „ 



Thus we obtain (2.18) and (2.19). This completes the proof of the lemma. 
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Lemma 10 There exists e small such that for \\(j>\\ Yl < £o an d t > 1 we have 
the following decay estimates 



sup 

0<s<t 



dip 
dx 



(s) 



< C 



L°°(x) 



d 2 ip 


(s) 


dxdx n 


d 2 ip 


(*) 


dxdxo 



ds < C- 



L°°(x) 



< c- 



t 2 



(2.23) 

(2.24) 
(2.25) 



Proof. Differentiating (2.13)-(2.14) with respect to x we get 



_9_ 

dx 



g (0 = v 2 (x {0 , |^ (0 = v 2 (x (0 , o [f J + ^ g («) ds] 



and thus 



< 



K + l) 

U\\ Yl 



N 



< 



< 



K + i) 

My, 



JV-1 



7o ^ (s)ds 



< +. 1 



L°°(a;) 



JV-1 



K + (e + 1) io 



dip 
dx 



(s) 



ds 



1 1 p 



L°°(i) 



Differentiating (2.12) with respect to x and using the estimate above, we obtain 



d<p 
dx 



(s) 



<c\\4>\\ 



L°°{x) 



1 



+ 



K + l) 

t 



JV-1 



JV-1 



* K+i)y 

i _l /■ 

* + k+i)7o 



9a; 



ds 



L°°(x) 



(is 



L~(x) 



< C — - + C \\4>\\ Y sup 

t ' 0<s<t 



d<p 
dx 



+ C\ 



log (i + 1) 



It then follows that, for t > 1 : 

5<p 



+ C\ 



(s) 

log (t + 1) 



sup 

I 0<s<t 



L°°(x) 



sup 

0<s<t 



<9x 



<C -^ + C||4 la sup 

L~(x) 1 ' 0<s<t 



dip 
dx 



(2.26) 
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By the assumption, if e is small enough, then we get 



sup 

0<s<t 



dp 
dx 



» 



< C 



Y,, 



L°°(x) 



and (2.23) follows. 

In order to derive (2.24) and (2.25), we compute -£^ a using (2.13), 

(2.14) and (2.15): 

(0 = v»* { x (o , o (0 JL* (o + ( x (0 , o J£_.y (0 
= ( x (0 , o [f / + / 1 («) «] [(i - 1 ) J + ^ ^ (,-) *] 

Jo 



dxdxQ 

1 d 2 G 
dxdxo 



(s) ds, 



(2.27) 



(o* + ij( { S5; ({) * (2 ' 28) 



Taking the norm IHI^oo^) of this equation, integrating the resulting formula 
with respect to s, using Lemma 9, (2.26), and the definition of \\4>\\ Yl > we 
obtain 



/' 

Jo 



dxdxr 



(s) 



ds 



L°°(x) 



< 



lo is 



d 2 G 



dxdxc 



(0 



< c- 



+ c 



t 

H\\ Yl 



t 



lo is 

I 



Jo 



d 2 G 



die 



<c^ + c,|,| k 



JV+1 



o 



dxdxo 



Yl.a 



(s) 



dxdx( 

1 d 2 <p 



(0 



L°=(i) 



dxdxo 
ds. 



dxdx 
(s) 



(s) 



ds 



L°°(x) 



ds 



L°°(x) 



o K + i) 



JV 



L°°(s) 



Thus if £o is small enough, we get 

r * aV 



/ 

Jo 



dxdxo 



(s) 



L°°(x) 



ds < 

t 
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and (2.24) follows. We now set s = t in (2.28) and use (2.27) to obtain 



d 2 <p 



dxdx 



(*) 



< 



L°°(x) 



d 2 G 



dxdxo 



(0 



dt 



L°°{x) 



< 



c\\n 



t 2 



C\ 



2 



N+l 



C\\4>\\ Yla C\\4>\\y la , 

" — + t 2 ~ C ' 

and the proof of Lemma 10 is complete. ■ 



* Jo a + ir -A> 



dxdxQ 



t 2 



In order to complete the proof of Proposition 3, it only remains to obtain es- 
timates for the Holder scminorms of ^ , , a d ' Jf . These bounds are obtained 

ox ' oxo 1 axoxo 

using ideas analogous to those used in the two previous lemmas. 

Lemma 11 There exists e small such that for t > 1 and \\<j>\\ Yl a < £ o, we have 
the following decay estimates: 



sup 

0<s<t 



dip 
dx 



(s) 



< C- 



0,a,(x) 



t^ c 



dp 
dx 

dip 
dx 

4 r 9V 



L 



dxdxo 

d 2 cp 



dxdxo 



(s) 
(*) 
(*) 
(*) 



0,a,(x) 



ds<C^ 
t a 



< c 



WW 



Y,. 



0,a,(x) 



t l+a ' 



(2.29) 

(2.30) 
(2.31) 



ds < C- 



0,a,(x) 



< C- 



t l+a 

Y 1>a 



0,a,(x) 



t 2+c 



(2.32) 



Proof. Using Lemma 8 and (2.23), we get 



0,a,(x) 



L°°(x) 



Jo d x 



(s) ds 



CUL, r£ 



ds 



< 



K + i) 
c 



f 

if Jo 



ds - 



L°°(x) 

cuw 



0,a,(x) 

f'+/ £ 

1 JO 



r 6 84>_ 

Jo dx 



L°o{x) 



Y 1 , a re 



0,a,(x) 



(£+l f +a t<* t 



< 



Yi,, 



(e + if h 



c r 



dx 



(s) 



ds + 



C\ 



lYi,, 



0,a,(x) 



(£ + l)^ -1 
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Differentiating (2.12) with respect to x, taking the Holder norm, and using the 
previous estimate, we get 



dx 



(s) 



< 



0,a,(x) 



< 



+ 



f 

J s 


1°® 


0,a,(x) 


C\\< 






t 1+a Js 










t 2+a Jo 


(e+lf- 1 



1 

t 

+ c 



d_ 

dx 



G(0 



d£ 



0,a,(x) 



C I 



+ 



sup 

0<s<i 



sup 

t 0<s<t 



(s) 



0,a,(x)Js (C + l) 

4 edt 



N 



r pa 

0,a,(x)Js (£ + 1 



N 



1 0<s<t 



90 
dx 



(s) 



where we have used that 



0,a,(x) 
L°°(x) 



0,a,(x) 



0,a,(x) 



as well 



< ^. We then deduce (2.29) 



as the fact that, due to (2.23), f I + §| 

provided eo is small enough. We now derive the Holder estimate of J^- (t). By 
interpolation, the Holder inequality, Lemma 9, and Lemma 10, we obtain the 
following two estimates 



/' 

Jo 



ox 



ds 



0,a,(x) 



< c 



f 

Jo 



1^0 {S) 



L°°(a;) 



d 2 <p 



< c 







dip 
dx a 



(s) 



dxodx 

\ 1— a 

da 1 



L°°(x) 



dxodx 



(s) 



ds 



L~(x) 



<C| 



,.1-q IMj 

Hii.a £a — fa 

< c 



0,a,(x) 



IT® 
ox 



L°°{x) 

Therefore, (2.30) and (2.31) follow. 



dxodx 



(*) 



< C- 



L°°{x) 



— a f2a 



< c- 



t 1+c 
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We now use Lemma 8, Lemma 9, Lemma 10, and (2.29)-(2.31) to get 
d 2 1 

J 0,a,(x) 



dxdx 

<c||v 3 </>(x(0,0 



lL°°(s) 



+ c||v 3 0(x (0,01 

+ ^[V 2 ^o, Q ,(,) 



(1 -I )/+ / 5S 



90 
9a; 



- 0,a,(a;) 



defy 



t 1 + J dx 



t 1 + Jo dx 



t 1 + Jo dx 



+ c||vV(*(0,0Loo (x) / 



a 

L«>(i) 

• ay 

dxdxo 



9x 



dxdxo 



is) 



ds 



J 0,a,(x) 



9x 



(5) 



0,a,(x) 



< 



< 



K + i) 
K + if 



V / 

Jo 

/ f 

Jo 



dxdxr 



€ r fl2 



9V 



dxdxo 



(s) 



(s) 



ds 



0,a,(x) 



ds + 



cu\\ 



Vi,. 



0,a,(x) (^ + l) JV t 1 +' 



Similarly, we get from (2.12) as well as the estimate above 
9 V 



dxdxo 



is) 



(2.33) 



< 



if- 



dxdxo 



0,a,(x) 



G(0 



0,a,O) 



«+ 7 / c 



9 2 



dxdxo 



G(0 



< 



+ 



4 ^ C I^Hn,o 



AT 



t 2+a Jo (e + l) W * Vo (£ + 1) 



(£ + 1) 



0,a,(x) 



9V 
dxdxo 

£ r 

dxdxo 



is) 
is) 



ds 



0,a,(x) 



ds 



0,a,(x) 
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By integrating (2.33) from s = to s = t, we have 



/' 

Jo 



t r & 



fly 

dxdxo 



(s) 



ds < 



c 



0,a,(x) 



+ j^ 2 - log (t + 1 



-^\o g (t + i) + cu\\ Yia j o 

1 r d 2 <p 



d 2 v 

dxdxo 



c\\<t>\\ Y , 



dxdxo 



cis 



0,a,(x) 







dxdx 



o 



ds. 



If e is small enough, then we obtain 
'* r d 2 <p 



o 



dxdxa 



ds < 



C I 



0,a,(x) 



t l+c 



Now putting (2.34) into (2.33) with s = t yields 



dxdx 



(*) 



< 



CM 



n, 



0,a,(x) 



t 2+c 



and this completes the proof of the Lemma. 



0,a,(a;) 



(2.34) 



2.3.3 The proof of Proposition 4. 

Now we prove the decay estimates for the higher order derivatives of <p. We 
prove Proposition 4 by induction on I. The induction hypotheses consist of the 
following estimates, for < m < £, 



sup 

0<s<t 



d m ip 



dx n 







dx m dxo 
d m+1 ip 



(s) 



C I 



< 



L°°0) 



sup 

0<s<t 



9"V 



ds < 



C I 



L~(x) 







t r ^m+l^, 

dx m dxo 



C I 



< 



0,a,(x) 



(is < 



C I 



(2.35) 



0,a,(x) 



dx m dxo 



(t) 



< 



(2.36) 



t m+l 



dx m dxo 



(*) 



0,a,(x) 



C\\4>\\ Y 

< "V"'° . (2.37) 



Estimates (2.35)-(2.37) have been already proved for in = 0,1 (cf. Proposition 
3). We begin with the estimates of G = G v in terms of ip : 

Lemma 12 Let i > 2 be an integer. Assume the induction hypotheses (2.35)- 
(2.37). There exists e such that for t > 1 and \\<f>\\ Yl , a < e Cb we have the 
following 



dx e 



(0 



< 



CU\\y tia 

Y I 



CI 



L~ (x) ^ + 1) 



+ 



(£ + i) 



aTTT SU P 

0<s<t 



d e <f 
dx e 



(s) 



, (2.38) 



L°°(a;) 
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— (0 



cH\\ Yta cu\\ Yta rav ■ 



— (f) 



< 



dx e dx 



\Y tlt 



rt r 



,a,(x)"*^ (f + lf (F + Tf 

Proof. Taking of G yi 



dx e dx 



0,a,(x) 

(2.39) 
ds (2.40) 

>(*) 



0«> 



yields 

l<i<£. 
Ji + ...+j<=< 

l<m<i 



cis. 

0,a,(x) 

(2.41) 



(2.42) 



l<m<i 



where 

X = X(0=z + ^-^+ / 
* Jo 

and where A i j 1 ,,j i are suitable numerical coefficients. Using the induction hy- 
potheses, we bound each term with j m < I for all 1 < m < i on the right-hand 
side of the above identity as 

a* ^, vt v^ v a* cu\\ Yta cu\\ Yi j cu\\ Ye j cu\\- 

n? :o>x - v ~ x x LOO(x) - ( ^ + 1)^-1 £T^- 

r a s 

< 



^ llfflly^ 



<9< 



t,(x) 



0,a,(x) 



E 

l<m< 

< c ™~ 

(f + 1) JV "V+ Q 

Putting the above inequalities into (2.42) yields 

^ ( ) ^ ( ,)-^(C + if- 1+ K + i) 3 



- E 1*^(^.0 

G " 

< 



j^- . . x 

ax-'™ 



L°°(x) 



I,oo (a:) 

n 



g x j P 



L°°(x) 



inequalities into (2.42) yields 

cH\\ Yta cu\\ Yla /•« av 
cw\\ Yt , a , CM*.. _ 



(is 



1) (4 + 1) o<s<t 



aV , 



>(*) 
L~(x) 



21 



d*G 1 CU\\ Y t 



.... , ddk. ^pV r -J 



< 



^ + l) JV - 1 + (^ + l) JV - 1 o| U ,< t 
In a similar manner, we take gfij^ 
8 e+1 G 



),a,(x) 



.to? (s) . 



0,a,(x) 



l<i<£. 
ji + ...+j<=< 



of G to get 



dxi* X 
dx e dxo 



j m <t 



l<m<« l<m<i 

We use the induction hypotheses to bound all the terms with j m < i for all 
1 < m < i on the right-hand side of the above identity as 

dx l ' dx d3 



8x~^) 



C 



) ddh x dji X 

dx dxh dx^ 

dx dxh W 



L°°(x) JO 



< 



..\\ Yla cu\\ Yta c 



L^(x) 
L°°(x) 

c 



ds ■ 



dxi* 



L°°(x) 
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Qi QQh 

dx 1 ^ ' ^ dx^dxi 1 



X ■ 



X 



0,a,(x) 



< 



dx 

Qi 

dx 



-v0(x,o 



0,a,(x) 



dd jl (p 



dx a dx^ 



(s) 



ds 



L°°(x) 



L=o(i) \J0 



« r 



99 J v 



dx dx^ 



jt(s) 



(is 



0,a,(x) 



+ £ 



2<m<i 



9* 
9~? 



(/ 



dxodxh 



n 
n 



9^ 



9a;- 7 '' 



Q X 3 P 



X 



< 



+ 



L°°(x) ) 



L°°(x) _ 
L°°(;r) , 



0,0:, (a:), 



(£ + l) i+Ar - 1+a i a tJ'i ' V* " fr 7 '* 

cU\\ Yln cu 



V 2 



lis, 



i+JV-l 



< 



Wily,, 



#1 



#2 



i- 7 * 



(e + i) iV 

Thus we obtain 
9 £ G 



dx 1 

d l+1 G 
dx e dxQ 



(0 
(0 



< 



C||0|| 



C U\W t , a 



< 



0,a,(x) " ^(C + lf (C + 1) ,V ^0 

This completes the proof. ■ 
We now prove Proposition 4. 

Proof of Proposition 4. Taking of ip and using (2.12), we get 



K + i) 

c 

+ 



JV 



Jo 



9 £+1 



dx e dx 



(s) 



ds 



L~(x) 



dx e dxo 



(s) 



0,a,(x) 



ds. 



9V 



9x £ 



t at 



-*5® = - 1 / e— (0^ 



9 £ G 



9x £ 



will denote a generic norm (or seminorm) that will be assumed to take 



the specific values 



d i ip 
dx 1 



(0 



< 



/' 

J s 



\L°° ' I JO, a • 



Then 



dx 1 



(0 



I/' 

* Jo 



d l G 
dx 1 



(0 



d£ < 2 



/' 

Jo 



9 £ G 
9a/ 



(0 



n 

1 Py£m,l 



dip 

dxi- 



-X 



L°°(x) 



(2.43) 

where we used | < 1. Suppose that the induction hypothesis (2.35)-(2.36) is 
satisfied. By Lemma 12, we then have 



— (0 

9o^ lU 



-(^ir T{ ^ + o s ^ 



9V , x 



h 



(2.44) 
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where 

Putting (2.44) 



if || 

+ a, if 



into (2.43) yields 

<c\\4>\\ Ye:a (J* — ^ 

Thus if e is small enough, we obtain 



■]o,a 



(2.45) 



^(0 



ox l oxo 



dx 

Using Lemma 



< 



(C + l) / * °< s <* 
(2.9). In a similar 

dx e dx (q} 



1 /■* 



way, we have 



Ivi « {77+ SU P 



0<s<t 



12 yields 



dx l dx KZ> 
i as in (2.45). 

to s = 



< 



;ing (2.47) into (2.46) and integrating from 



ds}, (2.47) 



where 7 is as in (2.45). By putti 
i s = t, we have 



s = 



dx e dxo ^ ^ 



ds<C 



i-t i-t 



Y t,. 



's {t + if J x f f Wo 



dx e dxa 



ds)} 



<c\\<p\\ Ye ,J^ + 

Thus we obtain (2.10) provided eo is small. We 
and put s = t in (2.46) to get 



dx e dxQ 



Qt+l 



dx e dxo ^ 
then substitute 



9x 

da)}. 
(2.47) for 



(2.10) 



Therefore the 
2.4 



1 

< - 



ft 



Jo {s + if - v +1 ■ 

proof is complete. 



2.4 Estimating the potential in terms of the density p. 

The following is a standard regularity result for the Poisson equation. 

Lemma 13 (Elliptic regularity theory) 
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Proof. For any fixed t>0we define 



p{z,t) = {t + lf p{z{t + l),t) 



Notice that: 



J \p(z)\d N z = J (t + lf\p(z(t+l),t)\d N z = J \p(x)\d N x 



On the other hand 



Ell v ^(^)llL~ (R «) = (*+ 1 ) Ar E( t + 1 ) £ ll v ^(-' t )l 



sup 

zi,z 2 GR JV 



\V k p( Zl ,t)-V k p(z 2 ,t)\ 



\z\ - z 2 \ 



1=0 



sup 

x t ,x 2 &L N 



L°°(R N ) 



\W k xP (x 1 ,t)-V k x p(x 2 ,t)\ 

i \ot 
\X\ - x 2 \ 



Then: 



f \p(z)\d 3 z+J2\KH;t) 



V k p( Zl ,t)-\7 k z p(z 2 ,t) 

i \ a 
\zi - z 2 \ 



< C Mx k , a , 
(2.48) 



On the other hand, by assumption 

A x <j) = p 

We define 



4>(z) = (t+l) N - 2 cf>(z (t+l),t) 



Then: 



A z cj) = p 

We now claim that the following estimate holds 



fc+2 



+ sup 



V*+ 2 0(*i,t)-V*+ 2 (*,,*) 



< CJ 



k 

I \p(z)\d N z + J2\\vip(;t) 



\zi - z 2 



li~(R«) + SU P „ 

^ ' Zl ,z 2 m N 



(2.49) 

\y k p(z 1 ,t)^v k p(z 2 ,t)\ 

\zi - z 2 \ a 



Indeed, a standard interpolation argument yields 



I LP 



< 



L 1 



< J , 1 < p < oo 



Using then the Calderon-Zygmund inequality it follows that: 



LP 



<C\\p\\ LP < CJ , Kp<oo 
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Therefore, the Sobolev embedding theorem implies that 



< q < oo 



Interior estimates for the Poisson equation in Sobolev spaces give a uniform 
bound on the W k+2 ' q norm of the restriction of (f> to any unit ball and hence 
of the C a norm of this restriction. Using this estimate, (2.49) follows from the 
inequality 



< C 



k+2,a-B 1 (x a ) 



C<*(Bi (*<,)) 



\\P\\k,a;B 1 (x ) 



< CJ 



that is just a consequence of classical interior estimates for the Poisson equation 
(cf. [7]). Using the estimate (2.48) it can be concluded that 



k+2 

E|| v ^(-'*) 

e=i 



sup 

z lt z 2 es. N 



<c\ 



\zi - z 2 



\x h „ 



Using the definition of <j> as well as the definition of the norm ||-||y as in 
(2.2), it then follows that 

U\\ YKa <c\\p\\ Xha 

and this completes the proof of the lemma. ■ 

2.5 Conservation of the L 1 norm of /. 

The following result is standard in the theory of the Vlasov-Poisson equation. 
See for instance [8]. 

Lemma 14 Suppose that f (x,v,t) solves the problem (1.1), (1.2). Then: 

|/ (x,v,t)\dvdx = / / \f (x,v)\dvdx , t>0 (2.50) 



2.6 The proof of Theorem 1. 



We now prove our main Theorem 1. We use a continuation argument as well 
as the assumptions on / to derive estimates for p and to close the argument. 
More precisely, we will show that an estimate of the form 



L~(R«) 



k 

f \p (x, t)\d N x + (t + if J2(t + if || VV (•, t) 

^u^^N+k+a \V k p(x,t)-V k p(y,t)\ 

+ (t + l) sup 1 1 7a 1 < Me 

K ,yGR« \x-y\ 
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for < t < t* implies an estimate of the form 

k 

\p (x, t)\d N x + (t+ if £ (t + I) 1 || W (-, t) || ioo(RW) 



. / , . -i \N -\-k-\-a \\7 k p{x,t)-V k p(y,t)\ 

+ (f + l) sup J : -a L < Ce (2.51) 

x, v em" \x-y\ 

for < t < t* + S (t*) , where S (t*) > and C is independent of M if e is small 
enough. Together with a suitable local existence theorem which guarantees that 
an inequality of the form (2.51) holds on some time interval [0,t\] for t\ > 
the estimates which have been derived imply that (2.51) can be extended to 
the whole time interval < t < oo and this implies Theorem 1. Notice that 
this in particular yields a global existence theorem generalizing that of [3] for 
N = 3. The local existence theorem can be obtained by combining the estimates 
obtained in this paper with a contraction mapping argument in a straightforward 
way. 

Notice that the decay assumptions on /o have not been used until now and 
only the decay properties of the potential <f> have been used. We now use the 
decay properties of /o for the first time in the following proof. 

Proof. Suppose first that < t < 1. Let (xo,Vo) denote the starting point 
for the solution of the characteristic equations reaching the point (a;, v) at time 
t. More precisely, 

^r = V{sh ^- = ^(X(s),s), X(t) = x, V(t)=v (2.52) 
Notice that X (s) = X (s; x, v, t) , V (s) = V (s; x, v, t) and: 

xo = x (x, v,t) — X (0; x, v, t) , v — v (x, v, t) = V (0; x, v, t) . 

Then 

h| = \vo(x,v,t)\ > \v\-J \Vcj>(s)\ds> M-C|4 t|> > \v\-^, (2.53) 
if eo is small enough. Taking the derivative -j^j of the formula 

P(x,t) = J f (x ,v )dv 

yields, for < t < k, 



1^ djpm q^q^(-0,V ) ..... dv, 

h+-+ji=e, 

a<i<l, 0<m<i 

Notice that the derivatives d ^-, are bounded < t < I &s C (l + \\4>\\ Ye ^ 
as can be seen by differentiating the characteristic equations (2.52) with respect 



27 



to x and v. It is then straightforward to see that for < I < k, using (2.3) and 
(2.53) yields 



dx l 



(*) 



< C 



Loc( x) 



< 



< 



(l + ll^llYi., 

Ce (l + 
C f e (l + U\ 



d\fo 



Q x i~mQ v m 

dv 

(i + H) K 
o < t < 1 



(xo,vo) 



dv 



(2.54) 



dx l [t) 



< 



0,a,(x) 



< 



Ce (l+U\\ Y( a ) J 
Cie (l + ||. 



dv 



(i + Hr 

a) ' - 1 ~ 1 



(2.55) 



Next we treat the case t > 1. By taking ^ of p(x,t) in (2.8), we get, for 
< ^ < fc, 



dx e 



(x,t) 



ji+...+ji+p=« 

0<i<£ 



av 

du 1 



(s ,V(0)) 



<9x 



v(oy... 



d \ n d p 

5) "<°»^ 



dct 



where 



X — Xq 

V (0) = V (0, t, x, w (t, x, x )) = — - h (f (0; t, x, x ) . 



By Lemma 9, Lemma 10 and the assumption (2.3), it is easy to see that for 
< I < fc, 



d l p 
dx l 



(x,t) 



< 



< 



ji+-+ji+p=e 

0<i<£ 



pi + -+ji+P+N 



1 



Ce (l + W 



Ye,c 



t e+N 

Using Lemma 8 with £ = k yields 



dxo 



(1 + M) 



K 



< 



dv 1 

Ce (l 



(xo,V(0)) 



In 



dx (2.56) 



N 



, t > 1 



^(i) 



< c 



in,, 



0,Q,(a;) 



< 



iJi + ---+j»+P+JV+a 

c £o (i + M ntC 



/ 
/ 



57o 



(xo,V(0)) 



+ 



dv* 



{xq,-) 



dxo 



J 0,a,(u) 

(2.57) 



(top C ^ ( 1 + ll^lln, Q ) 

l„ l\-R" - fk+N+a 



(i + kol) 



* > i 
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On the other hand, combining (2.50) with the decay assumptions for f in 
(2.3), it follows that 

\p(x,t)\dx < Ce (2.58) 



It then follows from (2.54)-(2.58), as well as from Lemma 13 that 
\\p\\x k , a < Ce + Ce U\\ Yk a < Ce + Ce \\ P \\ Xk a . 
Choosing e small enough, Theorem 1 follows. ■ 

3 CONVERGENCE TO THE SELF-SIMILAR 
BEHAVIOUR. 

We define the following set of self-similar variables 

f(x,v,t) = 1 g(y,v,r), (3.1) 
(t + 1) 

<f>(x,v,t) = * 2 $(y,v,r), (3.2) 

where 

y=JT^Y ) > r = Iog(t+l) . (3.3) 
A straightforward computation yields the following transformed system 

g T + (v-y)- V v g + e-( N -V T V y $ ■ V v g = Ng, (3.4) 

= J g (y, v, r)dv = p (y, t) , (3.5) 
where g (x, v, 0) = 50 (x, v) = f {x, v) = f (x, v, 0) . 

{r k \V k o(v t) — V k d(v' t)\ 1 

/ |^,0|^ + £||V^(.,i)|| + Sup ' |/ .,a ^ '\ 0<«<1 , 



l $ H^, Q = ™P L II V $ (•'*)IIloc (rW) + sup J -p L 

*>o 1 „ y^'eR N \y y 1 



Notice that Lemma 13 is also valid in self-similar variables 
Lemma 15 (Elliptic regularity theory) 

wn Yk , a <c\\p\\ Xk , a 
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We reformulate Theorem 1 in self-similar variables 
Theorem 16 Suppose that g (y, v) satisfies the following estimates 



EE 

e=o m=0 



d e go 



Q x mQ v l-m 



< 



(l + \x\) K (l + \v\f ' 



(3.6) 



E su p 

m=oX ,x'em N 
k 

E su p 

m=0 \v'-v\<l 



g x mQy°k- m (x,v) dxmd l°k- m (%', v) 



I il a 

\x — x'\ 

(x, v) dx m dv k ~ 



< 



(i + H) 



K 



, < a < 1 , 



dx rn dv k 



(x,v>) 



< 



\v — V 



{l + \x\) K {l + \v\) K 



< a 



(3.7) 



where K > N and e is small enough. Then there exists a corresponding solution 
of the rescaled Vlasov-Poisson system with 



\\p\\x h , a ^ °kS0 



The main theorem that we prove in this Section is the following 

Theorem 17 Suppose that the assumptions of Theorem 16 are satisfied. Then, 
there exist 9oc (y, y ) e Cf oc , p x (y) £ Cf oc (~l L 1 (R N ) , (y) e Cf+ 1,/3 satis- 
fying 



-Nt 



9 {y, Vo, t) -> 5oc (y, 2/o) , in cf c 



p{y,r) -> Poo (y) , 



fe 

oc 



/or any < /3 < 1, as r— >oo. Moreover, we have 

IIPooHlI(rN) = I|5o|Il1(rw xR N) 

and we /lave i/ie following representation formulae for g^ 
9oo (y, yo) = 9o (yo, y + ^oo (0, y, y )) 

Poo (y) = J goo (y, yo) Joe (y, yo) dy 

A y $oo (y) = poo (y) 



(3.8) 



(3.9) 



as well as the limit formula, as r — > oo 
9 (y, v,t) -» y 5oo (yo, y) 5 (« - y) Joe (y, yo) dyo , in V (R N x R N ) (3.10) 
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where w^, (s; y, yo) is the solution of the following integral equation 

uoo (a; y, V°) = - I e" (Ar_2)? V,$ (y + (y - y) e"« + jf * e^" 1 ^ fa; y, y ) d V , {J d£, 

(3.11) 



and w/iere (y,yo) is given by 



Joe (y, yo) = lim 



det —/at + e' 



(T-,y,yo, 



Remark 18 Notice that (3.8) can be read in the original set of variables as 



< \ 1 - 
POM) ~ JnP° 

as t — > oo, uniformly on sets \x\ < Ct. 



Remark 19 The function (s; y, yo) is small for small densities. In partic- 
ular the representation formula (3.9) implies that the rescaled density function 
Poo (y) approaches the one associated to the free streaming case, defined in (1.3) 
if £o — > 0. Notice that this shows that the particular profile that describes the 
self-similar behaviour of the solutions depends very sensitively on the initial data 
go- This contrasts with the situation in the one- dimensional case where the lead- 
ing self-similar behaviour depends only on the mass of the initial distribution 
but it does not depend on any other information on the initial data go (cf. [4]). 
However, notice that it is not possible to obtain a closed form expression for 
goo in terms of g due to the fact that the function (s; y, yo) depends on the 
values of the function $ for any t G (0, oo) . 

In order to prove Theorem 17, we introduce some changes of variables anal- 
ogous to the ones used in the previous Section. 

Suppose that the characteristics starting at yo, vo reach the points y, v at 
time r and we regard v = w (y,yo,r) , v = w(y,y ,r) as functions of y, yo, 
and t and make the change of variables from v to yo to get 



dv = 



det 



dw 
dyo 



dyo 



p(y, T )= / g{y,v,r)dv= e NT g {yo,v )dv = / e Nr g (yo, w (y, yo, *)) 
JR3 J J 



det 



dv 
dyo 



dyo- 



The corresponding boundary value problem in the self similar variables (y, v, t) 
reads 

Y(r)=y, Y(0)=y . 
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In the absence of the field, we solve 



as as 
Y(t)=v, Y(0)=yo, 



which yields 



T >, x y-yoe . y-yoe~ 

V(s) = — — , Y(s) 



1 - e- 



1 - e" 



yo-y 

1 -e- T 



As in the previous section, we formulate the above as a perturbed problem from 
the free streaming one. 

V = V + w, Y = Y + C, 

§-»-c, £-.-<»-»>■*.•(*+ c). 

C(t)=c(o) = o. 

It is straightforward to see that 

u (s) = f* e -( N - 2 )«V,4> (Y (0 , j-^r | T e- (Af - 3)e (l - e"«) V,4> (Y (0 , d£, 

(3.12) 



-(«-£) 



where 



Jo 



Along the characteristics, we have 

dv _ ( 1 
<9y V 1 - e ~ T 



e~TI N + *± ( t ) 



The following result provides some decay estimates for the derivatives of u>, 
analogous to the ones derived in Lemma 9. 

Lemma 20 There exists £q small such that for any r > 1 and any function $ 
satisfying 



we have 




duj 
dyo 


f 

Jo 




du ' 


f 

Jo 


oyo 



ds <CW 



Vb,a - £ °' 

duo 



L oo {y) 



ds <C\\<f>\ 



0,a,(v) 



Yo.a. ' 



dyo 

doj 
.dyo 



(r) 
(r) 



L°°{y) 



<Ce- T \m Yoa , (3.13) 



< Ce' T m\ 



o,a,(v) 



Yi, 



(3.14) 
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Proof. The method of proof is similar to the one used in the proof of Lemma 
9. We take J- of (3.12) to get 



< v -V;.!.;0{^— 



-(£->?) 



-£ 



Taking the L°° (y) norm yields 

<CJ$lly 0Q [ T e-V f - 2 K{e-£ + e- T + [ e 

°'° Js Jo 

e" r r e -(^{e^ + e- T + / 
Jo Jo 



-tt-u) 



<9w 
<9?/o 



<9w 



^(y) 



+ <7||$| 



-tt-u) 



<9w 
dyo 



dn}d£ 



L~{y) 



duj 
dyo 



(V) 



dn}d£. 



L°°(y) 



Integrating the above inequality from s = to s = t and using e ^ ^ < 1, 
e~ T < for jy < £, £ < r and N > 3 yield 



<9w 

%o" (s) 



< C||$| 



a / 

Jo 



e" 2s (is+ ||$| 



L=°(y) ' Jo 

+ C||*|| V „ e- T r + C||*| 



dr\ 



L°°(y) 



div 
dyo 



d v \. 



Thus we have 



duj 
dyo 



(s) 



ds < 



L°°(y) 



(3.15) 



provided e is small enough. We now specialize to s = t and use (3.15) as well 
as the fact that TV > 3 to get 



duj 

dyo 



Jo Jo 



-(t-v) 



<C|l*lka e " T + C, H*lka e " 



div 
oyo 

duj 

dyo 



dn}d£ 



dt d v 



where we changed the order of integration. We thus obtain (3.13). Using (3.13) 
we deduce (3.14) in a similar way. Thus we complete the proof. ■ 

We also obtain the following estimates for the derivative of uo with respect 
to y. 

Lemma 21 There exists ea small such that for any r > 1 and any function $ 
satisfying 



1*1 



< So, 
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we have 



sup 

0<S<T 



duj 



<C\\$\\ Yo , sup 



<9w 
<9y 



L°°(jf) 0<s<r 

Proof. We take ^ of (3.12) to get 

1 - e"« 



^ (S) . 



<C||*| 



9j/ 



0,a,(y) 



-(£-u) 



<9w 
<9y 



(»?) dry}de 



+ X e " (JV " 3)€ ( 1 - e " £ ) v S*( y (o.o{Y 



1 - e"« /■« 



1 - e 

Since TV > 3, we have 



=-(«-'?) 



9w 



sup 

0<s<r 



duj 



<C||*| 



Loo (!/) 



-(iV-2)£ 



{1 + sup 

\0<s<r 



<9w 



+ C r ll $ lk«. e ~ 



<cr||$|| y +cr||$|| y 



sup 

yo< s <T 


<9w 


sup 

\0<S<T 


0y- {s) 



Jo 



"(«- 



"(«-'?) 



L°°(j/) / 



Thus we have 



sup 

0<S<T 



dy 



(s) 



<cm\ 



L°°(y) 



Yo.a ' 



(3.16) 



provided £q is small enough. In a similar manner, using (3.16) we obtain 



sup 

0<S<T 



doj 
dy 



(s) 



<C\\n Yo C e- {N - 2)i {l+( sup 
Jo 



du> 
dy 



(s) 



-(£- 



^dr)}d£ 



[1 + sup 

\ 0<S<T 



duj 

% (S) . 



0,a,(y)J JO 



<c \m\ Yo , a +c\\n Y0iO 



sup 

, 0<s<r 



9w 

#1/ 



0,a,(v)/ 



0,a,(j/), 



This yields the Holder estimate of ^ and completes the proof. ■ 

We present the following estimates for higher-order derivatives similar to 
Theorems in the previous section. 

Lemma 22 Let I > 1 be an integer. There exists e small such that for any 
t > 1 and any function $ satisfying 
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we have the following 
d e oj 



sup 

0<S<T 



dy 



, (-) 



< C 11*11^ o , sup 

L°°{y) ' 0<s<t 



d u 
dy e 



<C\\$\ 



0,a,(y) 



Y,.„ ' 



f 

JO 



dy e dy 



dy e dy 



is) 
(r) 



ds <C\\$\ 



L°°{y) 



Yt, 



f 

Jo 



d t+1 uj 
dy e dy Q 



(s) 



< Ce- T m\ 



L~(y) 



Y(„ 



d l+1 uj 

dy e dy 



(r) 



ds < 

0,a,(y) 

< c e - T m\ 



Ye.a 



0,a,(y) 



Yt ■ 



As a consequence of Theorem 16, we have, for any k > integer and < 
a < 1, 

||p|lc*:.«(R N ) ^ Ce , || ( i > || C fc+2,c«( R jv) < Ceq. (3-17) 

We now study the limit behaviour of the self-similar system (3.1)-(3.3). Indeed, 
the limit behaviour is asymptotically equivalent to the free streaming case. 

3.1 Proof of Theorem 17. 
Proof. We begin with 

cj (s; y, y , r) = - f eT^-^V^ (Y (£; y, y , r) , d^ (3.18) 

J s 

-r^— T C e- (Ar - 3) « (1 - e-f ) V„* (Y (£; y, y , r) , d£. 

i 6 ./n 



+ 



By using (3.17) and by the dominated convergence theorem, as r — > cxo, in C fc+1 , 

w (s; y, yo, t) ->■ w (s; y, y , oo) = (s, y, y ) . 
In particular, we have, as r — > oo, in C k+1 , 

vo (y, yo, t) = u (y, v (y, y , r)) = V + u) (0; y, y , r) -> y + 1^ (0, y, y ) . 
Thus, we have, as r — > oo, in Cj„ c , 

e~ NT g(y,yo,T) = g a (y a ,v a (y,y ,r)) -> 3o (y ,y + (0,y,y )) = g^ (y,yo) • 
Next, since 



. dv 
dyo 



-, — )I N + e T — it) 

l-e T J dy 



using Lemma 20-Lemma 22 yields, as r — > oo, in C fc , 

<9v 



Nt 



det 



<9yo 



Joo(y,yo) ^1 + O(e ). 
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We then apply the dominated convergence theorem to get, as t — > oo, in Ct 



loci 



dv 



P(y,r) = J .go(yo,y + ^(0;y,yo,T))e T det— dy 

-> J 5o (yo,y + ^oc (0, y, y )) Joo {y, yo) dy 
= Poo {y) ■ 

Using the elliptic regularity theory from the equation 

A y $ = p, 

there exists (y) e C^ 1,/3 , for any < (3 < 1, such that 

Ay^oo Poo • 

Taking the limit in (3.4) as t — * oo yields (3.8). Finally, notice that, given a 
test function tp (y, v) 



• {y, v, t) tp (y, v) dydv = J e Nr g {yo, v ) ip {y, v) dydv = 



R 3 xR 3 



dv 



e T 9o (yo,w (y, y , r)) tp (y, w (y, y , r)) det —dydy 

Oyo 

and taking the limit t^oowc obtain 

/ 9 (y, v, t) ip (y, v) dydv ->■ / g (yo, y + w» (0, y, y )) ip (y, y) Joo {y, yo) dy dy 

JR3XR 3 J 

which can be written in the sense of distributions as 

9(y,v,r)^ J goo(yo,y)S(v-y)J 00 (y,yo)dy as r ^ oo 

This yields (3.10), whence the proof is complete. 
Notice that in the limit case £q — > (3.10) reduces to 



9{y,v,r) 



x> (yo,y)dy 
J 9o {yo,y)dyo 



5(y-v) 
S (y — v) as t — > oo 
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